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■^C [ S- and P- wave DD scatterings are studied in a meson exchange model with the coupling con- 

stants obtained in the heavy quark effective theory. With the extracted P- wave phase shifts and 
the separable potential approximation, we include the DD rescattering effect and investigate the 
production process e+e~ — !> DD. We find that it is difficult to explain the anomalous line shape 
ri ■ observed by the BES Collaboration with this mechanism. Combining our model calculation and the 

O ,■ experimental measurement, we estimate the upper limit of the nearly universal cutoff parameter to 

I ' be around 2 GeV. With this number, the upper limits of the binding energies of the S- wave DD 

Wh, and BB bound states are obtained. Assuming that the S- wave and P- wave interactions rely on 

the same cutoff, our study provides a way of extracting the information about S- wave molecular 
bound states from the P- wave meson pair production. 
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CSl ■ I- INTRODUCTION 

f^ ' Recently observed charmoniuinlikc states, called X, Y or Z, have motivated heated discussions on their properties 

^D , (see Refs. [l|-l5| for the detailed review). All of them are above the DD threshold and most states are near thresholds 

of two mesons. Various interpretations have been proposed, such as tctraquark or S-wave molecular states, hybrid 

states, dynamically generated states, or mixing states of cc state and exotic components, while the possibility that 

they are just cc states has not been excluded yet [6|, |7|. Among these possibilities, the molecular interpretation is 

. ^ worth consideration, because the molecular states are expected to appear near the threshold. 

j^ . To understand whether the proposed molecules exist or not, one may study the bound state solution of the heavy 

meson- ant imeson systems. The meson exchange models are widel y u sed in describing interactions of two hadrons 

[8l-[l7| . Other approaches include the gluon exchange models [I^, [l^ the unitarized model [2^ , lattice QCD [2l[ , 



X 



and QCD sum rule formalism [22|. According to these theoretical calculations, it seems that the existence of heavy 
quark molecules is inevitable, especially for the isoscalar hidden bottom states, which awaits the future experimental 
confirmation. 

Among the hidden-charm meson-meson systems, the DD system is the simplest one. Since Z? is a pseudoscalar 
meson, the total angular momentum of the system is equal to the orbital momentum. The possible quantum numbers 
of the system are J^'-^ = 0++, 1 , 2++, etc. There is no mixing between different partial waves, unlike in the case 
of the deuteron or DD* , where for instance the S ^ D- wave mixings are important. For the other charm meson and 
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anticharm meson system, there may exist open charm decay channels which renders the analysis more complicated. 
So we focus on the DD system and discuss the isoscalar case within one meson exchange model intensively. 

The investigations in the literatures indicate that such a scalar DD bound/resonance state may exist. In fact, a 
DD bound state was obtained around 3.1 MeV in a quark-based model [ig]. A quasibound state was also found with 
the unitarized method |20l . |23| . The vector meson exchange results in a possible binding solution [9| . Our previous 
results cannot exclude its existence within the meson exchange framework [14| and chiral quark model [1^1, either. 
However, a very recent analysis indicates that the existence of a DD bound state is difficult to understand [2^. In 
this paper, we will reanalyze this issue. 

Despite these efforts, the interaction in the heavy meson systems is still poorly known. It is expected that the study 
of the scattering provides us with additional information besides solving the bound state problems. For example, 
the S- wave DK scattering lengths and phase shifts give us additional information about whether the molecular 
interpretation for Z?sj(2317) is reasonable or not [2^. The scattering of D and D* off the X(3872) may reflect the 
DD* interaction if X(3872) is a molecular state [2g|. To better understand whether the DD system may form a bound 
state, we will calculate the partial wave scattering phase shifts and the relevant cross sections. This is the first major 
part of the present study. 

In addition to the observation of these unexpected hidden-charm X, Y, Z mesons, the BES Collaboration recently 
announced an anomalous line shape of the e^e^ — >■ hadrons total cross sections in Ref. [23. The structure is slightly 
above the DD threshold and slightly lower than ■0(3770). A similar anomalous line shape was also observed in the DD 
production [28|. The di- resonance assumption is one possible choice to understand such a structure [29[. However, it 
is difficult to identify just from the cross section whether this structure is due to a bound state, a resonance or the 
final state interactions (FSI). More studies are required. 

The dominant decay mode of ^(3770) is the P- wave DD. The FSI effects [30, |3l| have been considered in 
understanding its large non-DD decay observed by BES [32, ISj] . Their results indicate that FSI has non- negligible 
contributions. If the DD interaction were really strong, the rescattering effect would also lead to the anomalous line 
shape in their production. In Ref. [sJl , part of the FSI effects in the DD production has been included. No anomalous 
line shape appears. Here, we will study the rescattering effects in the process e+e" — >■ DD based on the calculated 
phase shifts and the Yamaguchi separable potential approximation [35| . Therefore, the multiple scattering effects are 
included. This is the other major part of the present study. 

With a heavier meson mass, a smaller kinetic energy and nearly the same potentials according to the heavy quark 
symmetry, the bottom analogous systems are more interesting. We will extend our study to the isoscalar BB cases. 

We organize our paper as follows. In Sec. [Hi we present the relevant Lagrangian, the derived potentials, and the 
definition of the threshold parameters. The numerical results for the S- and P- wave DD systems are given in Sees. 
Illll and llVl respectively. In Sec. |Vl we consider the rescattering effect in e+e^ — > DD. In Sec. I VI I the results for the 
BB cases are presented. The last section is our discussion. 

II. THE EFFECTIVE POTENTIAL AND THRESHOLD PARAMETERS 

One pion exchange between D and D is forbidden because of the parity conservation. What we need are the 
couplings of the D meson with the light scalar and vector mesons. The relevant effective Lagrangian in the heavy 
quark limit reads [l3, S, |33] 

C = gaTr[HaH] - il3vTr[Hv>' p^H], (1) 

where the field H denotes the degenerate (0^, 1^) doublet 

H = ^-^[P*^l, + ^Pl,]. (2) 

and u^ = (1, 0, 0, 0) denotes the velocity of the heavy mesons. The vector meson field has the form 

P, = ^^^P,. h=\ p -^ + ^ i^*o . (3) 

V K*- K*° 4> J 

For the coupling constant g„^ we use the value derived according to the chiral multiplets assumption [331, 9<y ^ 0.76. 
The values gv = it^p/ Itt = 5.8 and /Sy = 0.9 are obtained from the vector meson dominance [36l . |38|. 



Prom the Lagrangian, one may derive the DD scattering amplitude and then the effective potentials from the a, p 
and uj exchange, which are all Yukawa type 

^^('^) = -||:^"""^ 

Mr) - -3^e-"-, 

iDTrr 

Here we consider only potentials for the isoscalar DD system. For the isovector case, one changes the factor (-3) in 
Vp{r) to (+!)• So the interaction for the / = system is more attractive. 

In principle, one may solve the bound state and the scattering problem by inserting these potentials into the 
Schrodinger equation directly. However, we will see unphysical results appear because all the mesons are assumed to 
be pointlike particles. The situation is very similar to the TV TV case where the very short range interaction is unclear 
[39| . For the realistic system of composite particles, a form factor at each interacting vertex is necessary. We will use 
the monopole type form factor 

, ^ h? -m? 

n,) - -j^^, (5) 

where A '--^l GeV is the cutoff, rn is the exchanged meson mass, and q is its four-momentum. The improved potentials 
are [13 
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Here we use one cutoff to describe the system. The case A — ^ oo gives the former potentials. We will call this case 
the point particle limit. 

After one gets the phase shifts with these potentials, the threshold parameters can be derived through the definition 

limfc2^+icot(<5L) = — , (7) 

where oq (ai) denotes the scattering length (volume) for the S (P)- wave interaction. With this convention, a^ is 
negative if the interaction is repulsive. For a system with attractive interactions, oi, > if there is no bound state 
and ql < when one bound state appears. 

III. THE S WAVE DD SYSTEM 

We mainly explore whether there exists an S-wave bound state. From the previous studies [10, UM, [3, we have 
learned that the numerical results are very sensitive to the cutoff A. Since the contact interaction of the form ~ d{r) 
does not exist, it is instructive to study the case A — >■ oo first. We solve the bound state with the potentials in 
Eq. ((3|). The a meson is a broad scalar resonance with strong coupling to / = tttt S- wave scattering states. Its 
mass is not definite but around 400^600 MeV j40l - |42| . We choose two limit cases, rric ~ 400 MeV and 600 MeV. 
Other parameters are rup = 775.49 MeV, m^ = 782.65 MeV, and niD — 1867.23 MeV. If one ignores the vector 
meson potentials Vp and V(^ (we label this case NV), there are no binding solutions. After the inclusion of the vector 
exchange potentials (we label this case VC), we get a very deep bound state with the binding energy around 980 
MeV, which indicates this is not a physical solution. A smaller cutoff suppresses the contributions from the short 
range vector meson interactions. A number around 1 GeV should be reasonable. In the following calculation, we take 
A =0.8, 1.0, 1.2, 1.5, and 2.0 GeV and compare the results (note the lower bound should satisfy A > ruu:)- 

We present the calculated S- wave phase shifts of the elastic DD scattering in Pig. [TJ In these diagrams, E is the 
energy of the D meson in the center of mass frame. According to the Levinson theorem, the phase shift approaches 
180° when _E — )■ if a bound state exists. The results in the figure indicate that a DD bound state is possible if the 
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FIG, 1: The phase shifts for the S wave DD scattering with various parameters. The upper (lower) two diagrams correspond 
to the cases without (with) vector meson exchange contributions. The left (right) two diagrams are obtained with 771^=600 
(400) MeV, The cutoff A is in units of GeV. 

TABLE I: S- wave DD scattering lengths in units of fm, NV (VC) indicates the contributions from vector mesons are omitted 
(included). The number of * in the table indicates existence of a bound state. The binding energies are given in Table Ull 
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short range attraction is strong, e.g. A > 1.5 GeV. If the cutoff is around 1.2 GeV or less, the S- wave DD bound 
state does not exist. 

With the obtained phase shifts, it is not difficult to derive the S- wave DD scattering lengths. We list them in Table 
m The negative scattering lengths suggest that one bound state exists, which can be understood from the diagrams 
(c) and (d) in Fig. [TJ The results for the case A — ^ oo are also given in the table. The phase shift goes up from 180° 
in the case of VG, so the derived scattering length is positive. 

We also revisit the bound state problem which was studied in Ref. [ij]. The binding energies are summarized in 
Table U 

It is straightforward to get the S- wave total cross sections from the phase shifts. We present them in Fig. [2l One 
may check the values at the threshold with the formula (Jl=o = iiraQ. This explains why the cross section at threshold 
with A = 2.0 GeV is smaller than that with A = 1.2 (or 1.5)_GeV in the case of VG. 

From the above results, one concludes that an S- wave DD bound state may exist when the short range attraction 



TABLE II: The binding energies for the different cases in units of MeV. The three values for the S- 
A — !> cxD correspond to the ground state, the first and second radially excited states, respectively. 
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FIG. 2: The S- wave total cross sections for the DD scattering with various parameters. The upper (lower) two diagrams 
correspond to the cases without (with) vector meson exchange contributions. The left (right) two diagrams are obtained with 
m^=600 (400) MeV. The cutoff A is in units of GeV. 



is strong. Unfortunately the behavior of the short range interaction is not completely understood. Whether the vector 
meson exchange interaction is important, or equivalently whether A is large, needs further study. The determination 
of a reasonable range for this parameter is one major task in this framework. A — >■ oo is not a realistic case because 
the binding energy around 1 GeV is too large to be explained by the p and w meson exchanges. In fact, if this is the 
case, the DD bound state is so compact that it should be represented by a 4-quark system, and it is not consistent 
with the DD molecular state. In the cases of a finite cutoff, a value larger than 2.0 GeV will lead to the binding 



energy more than 30 MeV while that around 1.2 GeV does not result in a binding solution. However, it is difficult to 
identify the reasonable range without further information. The extraction of the scattering length from lattice QCD 
simulations or experimental measurements will be helpful. One will see that the P- wave BB production is another 
observable to constrain the range of A. We will come back to this point later. 



IV. THE P-WAVE DD SYSTEM 
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FIG. 3: The phase shifts for the P- wave DD scattering with various parameters. The upper (lower) two diagrams correspond 
to the cases without (with) vector meson exchange contributions. The left (right) two diagrams are obtained with 771^=600 
(400) MeV. The cutoff A is in units of GeV. 



TABLE III; The P- wave DD scattering volumes in units of fni 
omitted (included). 
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The P-wavc centrifugal barrier makes the interaction weaker than that in the S-wave case. Actually, we do not 
find any binding solutions even in the point particle limit (A — >■ oo). In Fig. |31 we show the phase shifts derived 
with different parameters. The diagrams indicate that the attraction is not strong enough to form a bound state 
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FIG. 4: The P- wave total cross sections for the DD scattering with various parameters. The upper (lower) two diagrams 
correspond to the cases without (with) vector meson exchange contributions. The left (right) two diagrams are obtained with 
m<^=600 (400) MeV. The cutoff A is in units of GeV. 



nor a resonance. Similar to the S- vifave case, one may get the scattering volumes from the phase shifts according to 
the definition Eq. ^. We present the numerical results in Table lllll The stronger the attraction is, the larger the 
scattering volume becomes. The obtained P- wave scattering cross sections are shown in Fig. 21 

It is very interesting to note that there is a resonancelikc structure or bump in the cross section, although no bound 
state or resonance pole (where the scattering phase shift crosses 90°) exists in this channel. The structure appears 
when the total energy of motion is around 40^150 MeV, depending on the parameters. The P- wave interaction 
probably has effects on the production of ■!/'(3770) or DD. We will explore this issue in the following section. 

To compare with the future experimental measurements of the scattering cross sections, one has to sum up different 
partial wave contributions. Here we have performed the calculation up to P wave since the higher partial waves yield 
smaller contributions. Because the S- wave cross section is much larger than the P- wave one, the total line shape of 
the cross sections atot = aL=o + <^l=i is similar to that of the S- wave case. That is, the resonancelike structure of 
the P- wave interaction does not appear in the total cross section. 



V. THE RESCATTERING EFFECT IN e+e -^ DD PRODUCTION 

A. Formulation 



The above P- wave structure motivates us to calculate the DD production by including their rescattering effect, 
which may be helpful to understand the anomalous line shapes observed by the BES Collaboration [23, ^3, ■ The 
schematic diagram for this effect in the process e+e^ — > DD is plotted in Fig. \5\ The DD pair comes mainly from 



■0(3770). In a previous work [3J], the production near the threshold has been studied by including the single loop 
contributions of the intermediate DD^ DD* + c.c, and D*D* . In this work, we would like to consider the multiple 
rescattering effects using a nonrelativistic method. 




FIG. 5: The rescattering effects in the DD production process. 



One may find the detailed procedure to derive the formula in the Appendix A. Here, we only present the resultant 
production cross section aprod = o'l + cr2 ■ The first part gives the cross section without FSI, 
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g5/2 



X |f.f. 



(8) 



where ae = 1/137 and f.f. indicates the contributions from intermediate vector resonances and background fields 
[3J, |43, 13- Since we here focus only on the illustration of FSI effect in the following analysis, we consider one 
resonance 0(3770) and take simply the form of f.f. from Ref. [43l |. 
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where the coupling g^oD is defined through {D[pi)D{p2)\i}}{p,\)) = -ig4,DD^^^^ ' {Pi - P2){'2^i^Y5^{p - Pi -P2), Qc 
is the electric charge of the charm quark, the decay constant /^ of -0(3770) is given by (O|c7pc|'0(A)) = /^m^e)i , 
Tn^p (Ft) is the mass (width) of -0(3770), is a relative phase, and F^j^{s) is an effective form factor describing the 
coupling of the virtual photon with the DD pair. Here, we assume 
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with Fq as an adjustable constant J43| . One derives g.i^jjQ from the branching ratio of the strong decay and /^ from 
the leptonic decay width Fee. Note one should consistently consider the rescattering effect in determining g^pDO from 
the 1/^(3770) decay. Now the decay width also has two parts F = Fi + F2. We have F2/F1 = cr2/o'i at the peak 

The second part reflects the rescattering effect 



U'K'^al\]2{ReUxImU)W-[{ReUf-{ImUf]W'\ ^ 



9 s^mu 



{Wf + {W^f 



(11) 



We have adopted the Yamaguchi separable potential approximation [35| in deriving this formula and A is the coupling 
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FIG. 6: The reproduced P- wave DD phase shifts 5l=i and cross sections (jl=i in the four cases, (a): 1/ai =700.0 fm~"^, (b); 
l/ai=58.0 fm~'^, (c): l/ai=520.0 fm"^, and (d): l/ai = 11.5 fm~"^, respectively. The sohd hue denotes the results calculated 
in the former section. 
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constant. This approximation has recently been used to study the structure of X(3872) [4^- In the above formula, 

ReU = -ml^r / dkj——-—— 5-t(fc)Fp(fc), (12) 



27r2 Jq (fc2 + 2M2)(fc2_a2) 

ImU = -mDj-^=t{a)Fp{a), (13) 



47r\/s 
6^ 



W^^ = 1-7ZJ^/ dk^,^-,[t{k)]\ (14) 



W^ - -^[i(«)]^ (15) 

1 R 

^('^^ = (fc2+ ^2)2+^2(^2 +^2)' (16) 

where a2 = m]j{y/s — 2to£)) and 7^ means the principal value integration. We will determine A, /3, and B by 
reproducing the phase shifts and the scattering cross section calculated in the previous section (see Appendix B). We 
have introduced a phenomenological form factor in the DD production vertex 

Mk) = ^, (17) 

where the cutoff Ap around 1 GeV reflects the loop contributions in considering rescattering effects. Here, k means 
the virtuality and the coupling between ■0(3770) and DD becomes small when the virtuality of the D meson is large. 
This form factor helps to derive a physically reasonable gj^oDi which will be obtained with the following expression: 

i3r(^(3770) ^ DD) ■ Tt = g^oD^ 24nm^ 



A 



9 mDm^, 



2{ReU X ImU)W'' - [{ReUf - {ImU)^]W' 
{wf + {W^f 

B. Results 



(18) 



We choose four cases of parameters to compare with each other: (a) nio- — 600 McV, A — 0.8 GcV, and NV; (b) 
m„ = 400 MeV, A = 2.0 GeV, and NV; (c) m„ = 600 MeV, A = 0.8 GeV, and VC; and (d) m^ = 400 MeV, A == 2.0 
GeV, and VC. The first case has the weakest attraction, while the fourth one has the strongest attraction. 

For the case (a), we have l/d = 700.0 fm-^. We found (;3,S)=(370,8.5), (380,5.2), (390,3.7), (400,2.7), and 
(410,2.2) [the unit of /3 is MeV] can all roughly reproduce the cross sections and the phase shifts. We take /3=390 
MeV, J5=3.7 as an example to illustrate the result. The reproduced phase shifts and cross sections are presented in 
Fig. lUa) and the calculated production cross sections are given in Fig. [71[a). We also put the BES data [28] in the 
diagram. When plotting the latter diagram, we have used: Ap=1.0 GcV, m^ = 3772.92 MeV [i^, Tt = 27.3 MeV 
[il, /^ = 100.4 MeV obtained from Fee = 0.265 keV [H], 5^,^^ = 12.6 from 5^(^^(3770) -^ DD) = 0.853 ^, 
Fq = 5.0 and (j> = tt/2. One finds that the FSI contribution is small in this case, and it certainly does not lead to an 
anomalous line shape. 

For the case (b), we have l/oi = 58.0 fm^^. The parameters may be (/3 (MeV),B)=(320,30), (330,10), (340,5.5), 
(350,3.5), (360,2.7), (370,2.0), or (380,1.6). As an example, we present the phase shifts and cross sections corresponding 
to /3=350 MeV, 5=3.5 in Fig. El^b) and the calculated production cross sections in Fig. W[b)- Now the coupling 
constant becomes g^^Q = 11.9, while the other parameters are unchanged. The attraction is stronger, but the line 
shape is still normal. 

For the case (c), 1/ai = 520.0 fm'^. The values (/3,B)=(370,50), (380,13), (390,7.0), (400,4.7), (410,3.5), (420,2.6), 
(430,2.1), or (450,1.4) are acceptable where the unit of /3 is MeV. With /3=400 MeV, 5=4.7, the fitted phase shifts 
and cross sections are plotted in Fig. IH^c). We present the calculated production cross sections in Fig. ^c). In this 
case, we have the same coupling constant g^j^Q = 12.6 as case (a). Because the contributions from the vector meson 
exchange interactions are small, the line shape is also similar. 

For the case (d), with 1/ai = 11.5 fm'^, one may use (/3,B)=(550,60), (570,12), (590,5.8), (600,4.5), (650,1.9), 
(700,1.0), or (800,0.2) to reproduce the phase shifts and the cross sections. Figure [Bljd) shows an illustration with 
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FIG. 7: The obtained DD production cross sections correspond to the four cases (a): 1/ai =700.0 fm^'^, (b): l/ai=58.0 fm^'', 
(c): l/ai=520.0 fm~'^, and (d): l/ai = 11.5 fm^'', respectively. We get the dash hues by ignoring the rescattering part (72 of 
the production cross section rj^rod = cri -f 0-2. The experimental data are taken from Ref. [231 . 



/3 = 650 MeV and B=1.9. We plot the corresponding production cross section in the last diagram of Fig. [T] Now two 
parameters have been changed: g^Df) = 5.6 and Fq = 2.0. For the other combinations of /? and S, the magnitude 
changes a little but it can be adjusted to the experimental data by varying Fq and the phase angle (/). Although the 
final state attraction is strong, the anomalous line shape does not appear even in this case. 

The results of these four cases (Fig. [7]) tell us that the DD rescattering effects cannot change the line shape of 
the production cross section. The reason is that the P- wave interaction is still not attractive enough. To see a cross 
section of a stronger attraction, let us arbitrarily consider an extreme case: 1/ai = 0.5 fm^'^. /3=500.0 MeV, and 
B=0.0, where no bound state or resonance pole is formed, though a sharp rise of the P-wavc phase shift is observed 
just above the threshold. The corresponding phase shifts and the resulting DD production cross sections are plotted 
in Fig. [51 We have used g^i^jj ~ 4.6 and Fq = 1.5. One finds that the anomalous line shape appears now. Also, the 
peak around 3770 MeV is shifted to a little lower position. The fact that the FSI may lower the mass of a bound state 
or a resonance reflects the couple channel effects. Because the scattering volume is much larger than the maximum 
number 0.1 fm'^ in Table [TTTl this case may not be realistic. 

From the above results, we conclude that although the strong FSI may lead to anomalous line shapes in the 
production processes in the unrealistic case, one cannot interpret those observed by the BES Collaboration with this 
mechanism. It is worthwhile to study further whether the anomalous line shapes are due to other nearby resonances, 
new resonances, or channel coupling effects. 
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FIG. 8: The P- wave DD phase shifts (a) and their production cross section (b) in an extreme case with l/ai=0.5 fm""^, /3=500 
MeV and B=0. We get the dash hne in the right diagram by ignoring the rescattering part 02 of the production cross section 
o"piod = cTi + o"2. The experimental data are taken from Ref. [28 |. 

VI. THE BB SYSTEM 

All the previous studies may be extended to the BB case naturally. Now T(4S') is slightly above the threshold. 
The features, in principle, should be similar to those of the DD case. But we will see more interesting results. We will 
present the figures only if the line shape is strange. When performing the numerical evaluation, the new parameters 
we have to use are the meson masses [H] m^ = 5279.34_MeV, mx = 10579.4 MeV, the width [4^ Ft = 20.5 MeV 
and Tee = 0.272 keV, the branching ratio Br(T(4S') -^ BB) = 0.96, the derived decay constant /t = 340.7 MeV, and 
the electric charge of the heavy quark Qi, = —1/3. The derived coupling constant without FSI is gxBB = 23.9. We 
adopt a little larger cutoff Ap=l.l GcV for this case. 



TABLE IV: The S- wave BB scattering lengths in unit of fm. NV (VC) indicates the contributions from vector mesons are 
omitted (included). The number of * in the table indicates that of the binding sohitions. The binding energies are given in 
Table nil 



rria 


Vector meson 








A (GeV) 






(MeV) 


exchange 


0.8 


1.0 


1.2 


1.5 


2.0 


CXD 


600 


NV 


0.026 


0.058 


0.083 


0.11 


0.13 


0.17 


400 


NV 


0.19 


0.26 


0.30 


0.34 


0.38 


0.44 


600 


VC 


0.039 


2.74 


-1.32(*) 


-0.46(*) 


1.08(*) 


-3.02(***) 


400 


VC 


0.21 


7.47 


-1.22(*) 


-0.38(*) 


1.76(*) 


-2.56(***) 



The S- wave BB bound state is more likely to exist than the DD one. Our calculated phase shifts do not exceed 
16° when one considers only the scalar meson exchange contribution, which indicates there is no bound state. After 
one includes the vector meson contributions, the phase shift starts from 180° with a cutoff A=1.2 GeV. It can go up 
to 210° from 180° with a stronger attraction A = 2.0 GeV. We present the derived scattering lengths in Table HVl and 
the binding solutions in Table [ill In the point particle limit, the observation that three solutions exist and that the 
binding energies are large indicates this case is not physical once again. A finite reasonable cutoff should be slightly 
larger than that in the DD case because of the heavy quark symmetry and the smaller size of the B meson. As for 
A=2.0 GeV and VC, the scattering phase shift crosses 6l=o = tt at a finite E and therefore the S-wave cross section 
becomes zero at that point. 

For the P- wave interactions, we show the phase shifts with various parameters in Fig. [9] It can reach 180° with 
a strong attraction, e.g. A =2.0 GeV, nicr = 400 MeV and VG. Note the line for A = 2.0 GeV in the diagram in Fig. 
9(c) is a resonance while that in the diagram in Fig. 9(d) is a bound state. Therefore, the results are more interesting 
than that in the DD case. We give the derived scattering volumes in Table |V] and the binding solutions in Table |lTl 
In the point particle limit, the positive scattering volume and the existence of one binding solution come from the 
observation that the phase shift goes up to ~ 245° from 180°. The large binding energies in this limit again require 
a finite cutoff. If the cutoff less than 1.5 GcV is reasonable, then no P- wave binding solutions exist. We present the 
cross sections for different parameters in Fig. 1101 Comparing this figure with Fig. 21 one observes the resonancelike 
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FIG. 9: The phase shifts for the P- wave BB scattering with various parameters. The upper (lower) two diagrams correspond 
to the cases without (with) vector meson exchange contributions. The left (right) two diagrams are obtained with 7TIct=600 
(400) MeV. The cutoff A is in units of GeV. 

TABLE V: The P- wave BB scattering volumes in unit of fm''. NV (VC) indicates the contributions from vector mesons are 
omitted (included). The number of * in the table indicates that of the binding solutions. The binding energies are given in 
Table [n 
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structure is more evident. It appears at a smaller energy of motion. Note that the cross section for the case A — 2.0 
GeV (c) also vanishes in the limit S = 0, which is understood with Eq. (jB7|) and the nonzero l/oi. 

To get the total cross section, we have to sum up all the partial wave contributions. In the DD case, the estimated 
maximum partial wave was Imax ~ ^'"^ ~ 1.87, so we calculated up to P wave. However, in the BB case, one 
should consider higher partial wave contributions (up to F wave) since Imax ~ 3.1. For the D- wave scattering, the line 
shapes look like those in the P- wave DD case, but the phase shift can go up to 100° with mcr=400 MeV, A — > oo, and 
VC. For a finite cutoff A < 2.0 GeV, a resonancelike structure also appears in the cross sections, although the phase 
shift does not exceed 90°. The phase shifts for the F- wave scattering do not exceed 40° even in the point particle 
limit. After summing up these four partial wave contributions, we get the total cross sections shown in Fig. [TTJ An 
interesting structure is there if one does not ignore the vector meson exchange contributions. The bump structure for 
the case A = 1.5 GeV comes mainly from the P- wave scattering. 
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FIG. 10: The P- wave total cross sections for the BB scattering with various parameters. The upper (lower) two diagrams 
correspond to the cases without (with) vector meson exchange contributions. The left (right) two diagrams are obtained with 
m<T=600 (400) MeV. The cutoff A is in units of GeV. 



Now we move on to the final state rescattering effects in the process e+e^ -> BB. One may still choose the four 
cases in the DD system. However, there are no strange line shapes in the cases (a), (b), and (c). Since both a 
resonance and a bound state are possible in the P- wave scattering with a finite cutoff, here we consider the case (d): 
rricr = 400 McV, A = 2.0 GeV, and VC; case (e): m^ = 600 MeV, A = 2.0 GeV, and VC; and an extreme case. 

For the case (d), we have 1/ai = —0.18 fm~^. Wc illustrate the results with /3 = 1200 MeV and B = 0.7. Figure 
[T^ shows the reproduced phase shifts and cross sections while Fig. [T51 shows the calculated production cross sections. 
We have adopted the coupling gxBB = 3-5 and the DD parameters Fq = 2.0 and (p = tt/2. The anomalous line shape 
reflects the P- wave BB bound state. For the case (e), we plot, with 1/ai = 0.22 fm^'^, gxBB = 3-0, (3 = 1000 MeV, 
and B = 3.0, the reproduced phase shifts and scattering cross sections in Fig. [Inland the corresponding production 
cross sections in Fig. [T31 One should note the production cross section vanishes at the threshold. So we have two 
peaks in the second diagram of Fig. 1131 The sharp one near the threshold is due to a P- wave resonance. Similar 
to the DD case, we also consider an extreme case: 1/ai = 0.5 fm~'^, /3 = 500 MeV and B = 0. Figure [14] displays 
the derived phase shifts and the obtained production cross sections where Qxbb = 11-3 has been used. Now the 
anomalous structure is just a bump. The peaks of T(45) in these cases are shifted to a lower position because of the 
strong final state interactions. 

The big difference between the two diagrams in Fig. [T2] results from the fact that the coupling constant in the 
separable approximation is sensitive to the scattering volume. To see the behavior of the line shape in the separable 
approximation, let us go to Eq. (jB5p . From that equation, the coupling is stronger with a bigger value of the scattering 
volume. If the scattering volume is positive, this indicates that the stronger the BB interaction is, the more obvious 
the anomalous line shape is. On the other hand, if the scattering volume is negative, the weaker the BB interaction 
is, the more obvious the anomalous line shape is. Therefore, if there is a sharp P- wave resonance or a shallow P- 
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FIG. 11: The total cross sections for the BB scattering with various parameters. The upper (lower) two diagrams correspond 
to the cases without (with) vector meson exchange contributions. The left (right) two diagrams are obtained with 771^=600 
(400) MeV. The cutoff A is in units of GeV. 



wave bound state around the BB threshold, the structure may be observable in e 



+^- 



BB. The sensitivity to 



the scattering volume explains why the two diagrams in Fig. [13] look so different. Since the value is positive for a 
resonance, the peak is clearer. In one word, the scattering volume mainly affects the line shape of the production 
cross section. The other two parameters /3 and B mainly control the magnitude of the cross section. But then one 
may recover roughly the same cross section by adjusting Fq and </>. 

In Ref. [43], the BaBar Collaboration measured Rbis) = cti,{s) / a fj^fj^{s) . One may derive a cross section around 1.2 
nb near T(4iS') while our result is around 2 nb. The reason may be that we considered only one resonance contribution. 
The interference between nearby resonances may reduce this number j44| . The normal line shape of the BaBar data 
tells us that the above three cases are not realistic, which means the P- wave BB interaction is not so strong. From 
this observation, we may get an upper limit of the cutoff, A < 2.0 GeV. In fact, investigation in detail reduces the 
upper limit to 1.7 GeV where there is a shallow P- wave resonance. 



VII. DISCUSSIONS 



From the calculated partial wave phase shifts, we know that neither the DD nor the BB bound state exists if 
one does not consider the vector meson exchange interaction. One gets the same conclusion from solving the bound 
state problem. Furthermore, in this case, the interactions for the isoscalar and isovector systems are the same, and 
therefore are not reasonable. The realistic interactions in the meson exchange models should include the vector meson 
contributions, and then the isoscalar interaction becomes more attractive. Our following discussions concentrate on 
the case 7 = and VC. 
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FIG. 12: The first (last) two diagrams show the reproduced P- wave BB phase shifts and scattering cross sections corresponding 
to the case (d) [(e)] with l/ai=-0.18 (0.22) fm"^, respectively. 
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FIG. 13: The obtained BB production cross sections correspond to the case (d) l/ai=-0.18 fm ^ and (e) l/ai=0.22 fm ^, 
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obtained BB phase shifts and production cross sections for the case l/ai=0.5 fm ^. We get the dash line in the 
by ignoring the rescattering part a2 of the production cross section aprod = a\-\- on. 



First, we focus on the cutoff A in our model study of DD and SB, which is to be almost universal for the bound 
state problem, the scattering problem, and the production problem. The sensitivity to the cutoff requires a reasonable 
range. This parameter for BB should be a little larger than that for DD and thus the upper limit from the former 
case applies to the latter case. Here we get an upper limit A < 1.7 MeV from the BB production cross section and 
the P- wave interaction. We may constrain the cutoff from the P- wave BB process because the possible resonance 
or bound state in the rescattering mechanism changes the line shape of the cross section. Since neither a P- wave 
resonance nor a bound state exists, one cannot get the constraint from the DD production. 

With the upper limit 1.7 GeV of the cutoff, we cannot exclude the possible DD S-wave bound state, but the binding 
energy is less than 10 MeV. There does not exist a P- wave BB bound state or a resonance, while the binding energy 
of the possible S- wave molecule should be below 100 MeV. We have used the vector meson coupling constants gv and 
^v derived from the vector meson dominance. The recent calculation with light cone QCD sum rule method gives a 
smaller gy ■ j3v [3- I* will lead to shallower meson- meson bound states. 

In our calculation, we have used a monopole type form factor in the meson exchange potential. One may alternatively 
choose a dipole type form factor. The potential in the latter case is weaker than that in the former one if the same 
cutoff is adopted. Provided one obtains the same scattering volume with Ai in the monopole case and A2 in the 
dipole case, A2 is larger than Ai. The corresponding S- wave binding energies have the relation \Ei\ > \E2\, but the 
difference is small. From the normal line shape of the BB production cross section and the P- wave interaction, we 
find that the cutoff now should be smaller than 2.4 GcV. The corresponding binding energies for the possible S- wave 
molecules change a little, but we still have the results: Ejjjj < 10 MeV and E^g < 100 MeV. That is, the usage of a 
different form factor has small effect on the constrain for the S- wave binding energies. 

If the DD S- wave bound state exists, its strong decay channel would be mainly rjcrj and Xcotttt. Other channels 
such as rjcTTTTTT and J /ipTTiriT are suppressed. This state may be produced in B decay, 77 fusion and pp collision [49|. 
Future precision analysis of the process e^e^ — > J/tp + scalar states may also offer a chance to find it. On the other 
hand, since the BB bound state is close to the thresholds of T{lS)(t> and Xho(2i')7r7r, its dominant strong decay 
channels would be T[lS)uj and Xbo{^P)'^T^ ■ One has a chance to obtain this scalar state by pp collision. 

In our approach, the S- and P- wave DD interactions depend on the same cutoff. Once the P- wave production 
including the final state interactions gives a constraint on its reasonable range, the S- wave interactions are better 
understood. This approach may be applied to new exotic resonances which arc candidates of S- wave meson-antimeson 
molecules. For example, the P- wave D*D* production will be helpful to understand Z+(4051) [501 ^'^'^ the relevant 
bound state problem |5ll . |52| |. The cutoff A will not change significantly because of the heavy quark symmetry. 
Another interesting example is the DD* interaction. The investigation on their P- wave production may deepen our 
knowledge about X(3872) ^^. 

From the study of the P- wave BB production by including the rescattering effects, we have seen the result in the 
separable approximation is sensitive to both the amplitude and the sign of the scattering volume ai. If there is a 
shallow bound state or a sharp resonance, the line shape of the production cross section may reflect that structure, 
which is difficult to identify just from the scattering cross sections. The S- wave production processes should have a 
similar feature. Since the S- wave system has stronger attraction, that process is more interesting. Such a study is 
expected to be helpful to understand some of the newly observed near-threshold structures. 

In summary, we have explored the bound state and the scattering problem for the isoscalar DD system, and the 
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rcscattcring effects in the e+e" — > DD process. We have also considered the corresponding bottom cases. From 
the binding energies and the phase shifts, the S- wave DD bound state would exist if the vector meson exchange 
interaction plays a major role. From the line shape of the calculated DD production cross section, it is difficult to 
understand the BES observation by the DD rescattering effect. From the line shape of the BB production cross 
section, we have estimated the upper limit of the cutoff A in the coupling form factor (monopole type) to be 1.7 GeV. 
Assuming this is the case, we get an upper limit of the S- wave binding energy: 10 MeV for DD and 100 MeV for 
BB. The future measurement of the S- and P- wave phase shifts, scattering length or volume, and the production 
cross sections may provide more information about the near-threshold resonances. 
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Appendix A: Production cross section 



Here, we illustrate the procedure to derive the production cross section after considering the rescattering effects. 
We use the PDG state normalization {p\q) = {^n)^ S^ {p — q) [4^. In the center of mass (cm.) frame, the DD system 
with the cm. momentum p also has this normalization. Therefore, one has 
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Our basic formula is the differential cross section for a 2 ^- 2 production process [46 
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where M is the Lorentz-invariant scattering amplitude, Ff is the flux factor and p (Ej) is the 3-momentum (energy) 
of the final state meson. 

To consider FSI, we insert Eq. (|XT|) into |Mp: 
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Here, the operator Oi describes the production of DD and O2 describes the rescattering of the final states DD. pmi 
and pm2 are the momenta in the cm. frame. They have different angles but the same amplitude. One may get the 
total cross section 
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where we have ignored the mass of the electron and have used Ff = s/2. 



d^Pf 
(27r)3 



m2\Ol\f){f\02\mi) 



(A4) 



19 



We calculate the first part as follows 



8e 



l\Ol\q)(p\0^\l) = —[-Ak-pk-q+p-qs] X |f.f.|^ 



(A5) 



where k is the momentum of the initial electron in the cm. frame and f.f. has been given in Eq. 
For the second part in Eq. (|A4[) . wc have 
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with G{E) = [E - H + ie]-^. When there is no FSI, G{E) = G"{E) = [E - H^ + ie]-^ , H" = 2{mD 
In general, 
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where Vint is the potential. 

To get the analytical expression of the cross section, we adopt the Yamaguchi separable approximation [351 . For 
the P-wave DD interaction, one may write down as 



{p\Vint\q) = 9{p) ■ g{q), 9{p) = Pt{p)p. 
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We determine A, /3 and other parameters in the function t(p) through reproducing the calculated phase shifts and the 
scattering cross section. After some calculations, one finally gets 
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In this formula. 
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With the equation p,_^^.^ — r-h ~ ^''^^{^ ~ ^)j it is easy to get the real part and the imaginary part 



E-H+in E- 



ReW = -mo-^V 

27r^ 

2^,3 



„ ^^-fc^t^W^ 



ImW 



-mo- 



/3^ 
47r 



■[t{a)]\ 



where a^ = ra£,{E — 2m£i). 

By combining Eqs. (jA4[) . (jAsp . and (jA9[) . we obtain the resultant cross section cr = cti + (T2, 
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where Qfe = 1/137 and 
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If one defines W = 1 + rr^ReW, W = ^ImW, then 



64 TT^a^A j 2(i?eC/ x ImU)W - [{ReUf - {ImUf]W' 
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When there is no FSI, A = and CT2 vanishes. 

Appendix B: The parameters 

Before the numerical evaluation, we have to choose the form of the function t{k) and determine the relevant 
parameters. In Ref. |55| . the separable P- wave potential was given with t{k) = ^, ^_/"g^.^ . However, we find the 
following choice is better when reproducing the scattering cross sections: 

1 R 

^(^) = (F+^2)2+^2(p+^2)' (Bl) 

where i? is a dimensionless parameter. One notes the second term of g{k) does not have a good behavior at large 
k. This is easy to see after the Fourier transformation. Since we consider only low-energy interactions, this form is 
acceptable. 

From the scattering amplitude, one has 

,■1 .. ^ 127r 4TTF(k) ,„ , 
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where Si is the P- wave phase shift. According to the definition of the scattering volume ai, 
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So the coupling constant is 

127r 
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If one does not explicitly use the coupling constant in expressing the P- wave cross section, we have 
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By inserting the function t{k) into the above formulas, it is easy to get the exphcit expressions which we do not 
present here. The scattering volume has been derived in Sec lIVI Now the parameters we have to determine are /3 and 
B. They are extracted by reproducing the Si and aL=i- 
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